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NOTE ON STELLAR PERTURBATIONS OF NEARLY 
PARABOLIC ORBITS. 


By E. Optx. 
(Read by title October 12, 1932.) 


TuHE following considerations refer to the question of statistical 
stability of cometary or meteoric orbits in the solar system, and partly 
also to the stability of wide double stars as far as perturbations by 
passing stars are concerned. The probability of ejection of a member 
out of the solar system, and the probability of getting out of the 
region of observability of the terrestrial observer will be discussed. 
Certain simplifications have been introduced, and on this account the 
estimates are approximate, being supposed to give the order of 
magnitude only. 

An object (which we will call the meteor for the sake of simplicity) 
of very elongated orbit will spend most of its time near aphelion, and 
the stellar perturbations are there the strongest and the most frequent. 
For the order of magnitude it is therefore legitimate to consider only 
perturbations caused when the meteor is in aphelion. Assuming for 
the calculation of the frequency of perturbations that the meteor 
remains all the time in aphelion, we slightly overestimate the effect 
of perturbations, without much affecting the order of magnitude; on 
the other hand, the probable density of stars in space may be under- 
estimated, because the number of white dwarfs, which may be rela- 
tively massive, is not known and is not taken into account here. 
This underestimate may balance the overestimate mentioned above. 

The following considerations are valid when the effective duration 
of approach of star and meteor, defined later on, is small compared 
with the period of revolution of the meteor. This condition is easily 
fulfilled in all cases of practical interest, because of the high average 
velocity of the stars compared with the aphelion velocity of the meteor. 
Neglecting very close approaches of stars to the sun, the path of the 
star, on account of its high velocity, may be assumed to be rectilinear, 
with respect to the sun as well as to the meteor; the latter, during the 
short time of passage of the star, may be considered as practically at 
rest. 
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By reason of symmetry, after the passage is over, all the gravita- 
tional impulses parallel to the path of the star will be cancelled, and 
the resultant impulse contributed by the star to the meteor will be 
directed perpendicularly to the rectilinear path, toward the point of 
closest approach. Let us denote by & the distance of closest approach 
of star to meteor, expressed in astronomical units; by v the velocity 
of the star relative to the sun, expressed in units of the orbital velocity 
of the earth; by a the angle between the line joining star and meteor, 


and the direction of closest approach. Choosing as unit of time 2. 


years, and assuming the mass of the star equal to the sun’s mass 
(which may be near the average), the gravitational acceleration ex- 
cos” « 


erted by the star upon the meteor will be ; the component of 
this force parallel to the line of closest approach is 
_ cos? « 


F, = 
4 


and the total impulse contributed to the meteor along this line is 


Au = 2 F 
0 


the time ¢ being reckoned from the moment of closest approach. On 


the other hand, 
v= where 2x = étga, 
dt 
whence 
dt = mit 
v COS” & 
This gives 
9 
Au = — f cos ada 
fv Jo 
or 
Au = (1) 


The geometrical sum of the original velocity of the meteor and of 
Au gives the perturbed velocity; the latter determines all the elements 
of the changed orbit. 
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The effective duration of approach, 4, may be defined as the time 
necessary to produce the whole effect, Au, by the star placed at the 
point of closest approach; thus 


1 
ty Au, 
whence 
v 


One of the conditions of validity of (1) is that the period of revolution 
of the meteor, P, must be great as compared with f. For a very 
elongated orbit we may put 


in the adopted units of time, whence 


« 
ov2 
This ratio must be great—at least ten or more. The other condition 
for the practical validity of our considerations is that the eccentricity 
be great, that is, greater than 0.9. 
For the average velocity of a star relative to the sun, we assume 


v = 1 = 30 km/sec., which is approximately the observed average. 
In this special case our formulae are 


(3) 


9 


Au = — 
(1) 


and 


3/2 
= 1,1. (3’) 
0 


For a highly eccentric orbit the aphelion velocity is very small com- 
pared with the velocity of escape at aphelion distance; when con- 
sidering the probability of the ejection of the meteor out of the solar 
system, we may neglect altogether the aphelion velocity, assuming 
it equal to zero. The condition for ejection is, then, that the relative 
perturbation, or the geometrical difference of the perturbation vectors 
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Au, and Au, for meteor and sun, be equal or greater than the velocity 


of escape, given by 
2 
e= — 4 


in the units chosen here. The condition of ejection is 
Au = », (5) 


Au being the absolute value of the relative perturbation vector for 
meteor and sun. It is easy to show that when R is great, it is neces- 
sary for ejection that the inequality 


Se 4/2 6 
n= R (6) 


be fulfilled; this follows when for a given value of £, we consider the 
case of strongest effect, when the star passes through the line joining 
m and s between them, and when & = R—£é,. Here it is assumed 
that &,, < &. When &, > &, we obtain the same result by interchang- 
ing the places of meteor and sun. When é,, > R (é > R also), the 
absolute sum of the two vectors 


2 2 
Aum | + | Aus =—+—<-, 
| Adm | + | Ares | 


is smaller than 2, for R > 8; thus ejection cannot take place. From 
(6) we infer that for R exceeding about 100 astronomical units the 


ratio se is small, whence &,, is small compared with R, and thus small 


relative to & also, and in (5) we may neglect Au, altogether. Ejection 
may happen only when the passage is so close to one of. the two 
bodies that the relative perturbation of the other may be neglected. 
This gives the following simple form for the condition of ejection: 


= 22 4/2, or fm S VW2R; 


bn 
the effective radius of approach for ejection is thus 
(7) 
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The validity ratio (3’) in this case is 


3/2 
to V2R 


which is large for all possible values of R that may be of practical 
interest in our problem. 

Assuming an average of twenty-five stars of mass 1 per 1000 cubic 
parsecs, or 3.10~'8 per cubic astronomical unit, and taking into ac- 
count that the total target area for ejection is twice zé2 (because a 
close approach may happen with respect to the sun as well as with 
respect to the meteor, and because these two target areas do not 
overlap when & < 14R, or when R > 8), we find for the probability 
of ejection within At years the value 


3.10718 - - (8) 


the last factor, 2x, is the number of astronomical units traveled by 
the star in one year. After substituting (7), the differential expression 
(8S) leads to the following finite expression for the probable fraction, 
ps, of meteors not ejected, having an average aphelion distance R a. u., 
after an interval of time equal to ¢t years: 


= 4.10°6Rt (9) 
This formula provides us with a means of estimating the maximum 
possible limit of the solar system. Substituting for ¢t a value near 
3.10°, the probable age of the solar system, we find for the probable 
fraction of meteors retained since the origin 


p, = 10°R, 

The effective limit, corresponding to a fraction 1/e retained out of 
the original number may be put thus at R, = 10° a. u. = 5 parsecs. 
The effective volume of the solar system is about 600 cubic parsecs, 
or about fifteen times greater than the volume per average solar mass 
in space. 

This result holds for the assumption of a constant density of stars 
in space; the average density of stars along the path traveled by the 
sun since the origin of the solar system might have been different 
from the present density; there may be a variation with time also. 
That is, if the galactic system has been expanding appreciably, the 
average density of stars in space must have been greater in the past 


5 
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than it is now, and the effective radius, R,, may be assumed to have 
been somewhat smaller. Taking into account this circumstance, we 
may estimate the possible limit of the solar system at from 2.10° to 
10° astronomical units. The period of revolution for these extreme 
distances will then be equal to from 3.10’ to 4.108 years; thus, even 
the most remote present members of the solar system must have had 
time to complete several revolutions since the time of their origin. 

The quantitative conditions of ejection, derived here for nearly 
parabolic orbits, will not change very much, as far as the order of 
magnitude is concerned, for less eccentric orbits nor even for nearly 
circular ones. Thus, the effective limits set here may be regarded as 
of general validity. 

For a star of mass uv. times the solar mass, the exponent in formula 
(9) is to be divided by w; for the same age, t, the extreme limit of the 
system R, is to be multiplied by yu. In this form our considerations 
may be applied to the question of stability of wide double stars—the 
so-called pairs with common proper motion. For uw = 1 we find the 
time of retaining the fraction 1/e of companions at an aphelion dis- 
tance R as follows: 


t, years 3.10° 3.10!° 3.10" 3.10” 


For other values of the mass these intervals of time are to be multi- 
plied by u. It may be shown that the average velocity of ejection is 
such that within the time ¢, for R < 10°, the separation of the two 
former members of the system becomes much greater than fifty par- 
secs or 10’ a. u. so that in most cases an ejection results in practical 
separation; the former components will not be counted any more as 
members of the same system. 

In the preceding it was assumed that the orbit, or at least the 
aphelion distance, remains substantially unchanged until a close ap- 
proach throws the meteor out. This assumption is not exactly true. 
Less powerful, but more numerous perturbations will have the effect 
of increasing the major axis and the average aphelion distance, thus 
making the final ejection more easy. The rate of ejections is therefore 
somewhat underestimated above; the order of magnitude, however, is 
not affected. As to the effective limit of the solar system, the estimate 
need not be changed, because the shifting outward of R will increase 
the limit and counterbalance the effect of a greater number of ejec- 
tions. The conclusion is that within intervals comparable with the 


Bis 
hrs 
a 
= 


STELLAR PERTURBATIONS OF NEARLY PARABOLIC ORBITS 175 


age of the solar system, meteors or comets remain practically well 
connected with the sun even at distances as great as the average 
separation between two stars in space. However, from the stand- 
point of observability their detection becomes impossible when the 
perihelion distance exceeds 1 for meteors, and is a small multiple of 1 
for comets. The perturbations of the perihelion distance are thus the 
most important ones from the practical standpoint. The question of 
persistence of gaseous constituents in meteors, in connection with 
the appearance of comet tails and with the determination of the age 
of meteorites from the helium content, also depends chiefly upon the 
perihelion distance, which determines the maximum temperature at- 
tained by the meteor during its revolution around the sun. 

The law of conservation of angular velocity gives for a nearly 
parabolic orbit, when the perihelion velocity is assumed equal to the 


2 
velocity of escape, 


Qn 2 
= (10) 
here q is the perihelion distance, R the aphelion distance, v, the 
aphelion velocity (units same as before). Thus q is very sensitive 
to changes in 2, being proportional to the square of the latter. The 
calculation of the perturbational vector in the present case is more 
complicated than in the former, because we cannot neglect one of the 
two impulses, Au, and Au,. Without going into detail, we may 
state that the absolute value of the relative impulse meteor manus 
sun is given by | 


Au = 


q 


2R-|sin | 
Em V m2 + Rk? sin? y + 2KE,, (cos « cos A + cos 6 sin A) 


here a, B, y are the angles formed by the vector R with the a, y, 2 
axes of Cartesian codrdinates; the z axis is assumed to coincide with 
the path of the perturbing star, and the meteor is placed in the zy 
plane; 1 is the azimuthal angle of the line of closest approach to 
meteor, recorded from the x axis in the xy plane.' The independent 


1 The derivation of (11) is simple in Cartesian coérdinates, where 


and 


Au = VaAu;? + Au,’? + Au,’. 
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variables are é,, R, and the direction of R defined by a, 6, and y 
under the restriction cos? « + cos? 6 + cos? y = 1. From (11) it is 
possible to find the frequency of different values of Au. The exact 
procedure, however, is too complicated; in this preliminary discussion 
we will substitute in (11), for the trigonometric functions, their 
average values calculated on the assumption of a random orientation 
of the vector R. In this way we obtain a certain approximate expres- 
sion for the average value Au of the vector au. For the function (11) 


this procedure means underestimating Au. F urther, let Av, be the 
component of Au along 2%, so that 


Av, = Au cos ¢, (12) 


where ¢ is the angle between Au and x. Assuming a random frequency 


for this angle,| we have | = 4%. Also, in (11), | sin y | = 
sin? = cos = cosB = 0. This gives for the order of magni- 
tude of the average perturbation in », the expression 
R 

Vint + 
From (10) we find the expression for the changed perihelion distance 


_ BR (vq + Ary)? _ 


2 2 


(13) 


Ae, = + 


R? 
+ > (Ata)? + Ara 
and with the aid of (10) again 
(Ana) +R V 2q Ary. (14) 


Substituting for Av, its average value from (13), retaining the + sign, 
and putting ¢ = &,,/R, we get the average perturbation of q as follows: 


32 + 24) % 


The first member is always positive and gives the average systematical 
increase of q as the result of stellar perturbations; the second member 


' This assumption is not exactly fulfilled, because ¢ depends upon a, 8, y: 
the vector vq is at right angles to the vector R. 
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is symmetrically negative and positive and cancels out in the mean; 
this member, being proportional to Av, which was assumed to be 
the projection of a random vector Au on v%, may be used to estimate 
the dispersion of Ag = qi — q around its average value given by the 
first member of the formula. Neglecting the dispersion in Au, we get 
a minimum spread of Ag, according to the frequency function 


sin gdg = dt, (16) 
where 
{= Ag = cos ¢ 
2Aqo 


and Aq is the absolute value of the second member in (15). 
Now, the main spread in Au, according to (11), depends upon the 


factor sin y, for given £ and R, or given< = , except when ¢ = | which 


is a peculiar case. ‘To allow for this spread, we have to put, instead 
of ¢, 


t 
fh = X constant; 
sin ¥ 


the frequency function of y is sin y. The frequency-function in- 
fluenced by the spread in Au will be 


y=2/2 
f dt; X frequency-function of + 
y=0 


-f —* sim ydy X constant = dt X constant. 
0 sin ¥ 


Thus we see, that the peculiarity of the spread in Aw is such that it 
does not much influence the spread in ¢ which is determined by (16). 
According to this formula, all values of Ag from — 2Aqo to + 2Agqo are 
equally probable. 

Equation (3’), upon which the validity of our formulae depends, 
inay be represented as 


0 


For practically important values of ¢ it may be shown that 7 isa 
0 


great number for R exceeding 1000 a.u. The extreme limit of a rough 
applicability of our formulae may be estimated at R > 300 a. u. 
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The target area for an approach at a distance = <R may be put 
equal to 


Fre? (17) 


c=F.- x? 


The factor F is equal to 2 for ¢ S 0.5 and decreases slowly for 
greater values of <¢ toward 1; this factor depends upon the size of the 
overlapping fraction of the areas of two circles, each of radius &, 
centered at sun and meteor, the distance between the centers being 
thus R. The formula for F is 


where 


The probable number of passages closer than & = ¢f within the 
time ¢ years is given in the present case by 


N, = 3.10718. Fxe2R?- 2xt = 0, 60.10- (8’) 


Table 1 gives the results of a sample computation for g = 1. 
In this table < = e, where é is the distance of closest approach either 


to meteor, or to sun, and F# is the aphelion distance; q; is the changed 
perihelion distance; the average increase, g; — q, does not depend upon 
q, the initial perihelion distance, unless g is not small as compared 
with f. The effective limits of g may be somewhat surpassed in the 
actual case, but for practical purposes it is safe to assume that the 
maximum and minimum values of q, given in the table, symmetrical 
with respect to “gq; average,” represent the extreme limits of variation 
of g, and that the frequency function of q is constant with respect to 
mi, that is, that equidistant intervals of gq are equally probable, to 
their extreme limits. The negative values of “gq: minimum”’ are 
introduced as a matter of convention, to indicate that the possible 
values include gq, = 0. Thus, from the table we infer that fore = 0.125 
and less, the changed perihelion distance is always greater than g = 1. 


ee or 0.003 will have |q:| less 
36.1 — 0.9 


For < = 0.16, only the fraction 


| 
than 1; fore = 0.25 the same fraction is ane or ten per cent, and the 
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fraction of perturbed orbits which will approach the sun to a distance 


closer than |q:| = 0.25 a. u. is i = 2.5 per cent. From ¢ somewhat 


less than 0.20 to about 1.3, a certain fraction of the perturbed objects 
may hit the surface of the sun, when | q: | < 0.0045; for ¢ = 1.28 this 


ato : or one in 250. The 6th and 7th lines of the 


table refer to data connected with N, as defined by (8’); the probable 
fraction of objects which in the indicated interval of time have not 
undergone an approach, defined by the given value of ¢, or closer, is 


(9’) 


Thus, from the table we infer that at R = 10500 the fraction 1/c of 
all objects has not suffered from approaches closer than ¢ = 0.16 or 
1600 a. u. within 3000 million years; or that at R = 2000 a. u. the 
fraction ¢~?:© or about 8 per cent has not suffered changes from ap- 
proaches closer than < = 1.6 or 3200 a. u. The table, though con- 
structed for g = 1, may be easily adapted to any other value of g, 
when g remains small compared with R; in agreement with formula 
(15), we must only add g — 1 to “gq; average,” and change the range of 
the effective limits of qg; in the ratio ¥/q. All the other data of the 
table remain unchanged. Any changes in JN,, R or t in the last two 
lines of the table are to be made in agreement with formula (8’). 

The most remarkable fact revealed by Table I is the tendency of 
the perihelion distance to increase continually. After a sufficient 
interval of time the perihelia will shift out of the sphere with radius 
q = 1, and even out of the known planetary limits of the solar system. 
The systematic increase in g, depending upon the probability of an 
approach, increases with the aphelion distance. Thus, for F suffi- 
ciently great, we may expect that most of the meteors, or comets, have 
already gone out of the region of observability. This limit of R may 
be estimated from the table at 10000 to 20000 a.u., assuming an age 
of 3.10° years, and would correspond to periods of revolution from 
400,000 to one million years. It is not very probable that many of the 
observable meteors have periods of revolution greater than these 
figures, unless the total number of such meteors in the solar system is 
very great. In the latter case, because the spread in q: increases as 
/q, a small fraction of meteors possessing a large g may be changed 
by perturbations into those with a small value of qg; and such objects 
may thus enter observable regions of the solar system. A state of 
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statistical equilibrium will be installed when the number with small qg 
changed into big q; will be balanced by the number with big g changed 
into small g;. Because of the said systematical trend of the perturba- 
tions in q, the equilibrium will take place when the number with great 
q will be greater than the number with small g. It is easy to see that 
this means also a certain increase of concentration of meteoric matter 
in the nonobservable regions of the solar system, although the con- 
centration is only a relative one, being counter-balanced by dilution 
in a greater volume of space. If, however, among our observable 
objects there is a sensible proportion of meteors with R exceeding 
10000 a. u., it may be only the result of a very peculiar distribution of 
meteoric matter at the outskirts of the solar system—a kind of a 
meteoric cloud, or shell, of a density much higher than may be an- 
ticipated from the observed frequency of solar meteors caught on the 
surface of the terrestrial atmosphere. Unless this is the case, we may 
expect that the majority of solar meteors we observe comes from dis- 
tances less than 10000 a. u. 

This supposition means that an average solar meteor must have 
completed at least many thousand revolutions during its life time, 
being thus exposed for a total interval of a thousand years or more to 
comparatively high temperatures during perihelion passages. This 
deduction may have some bearing on the determination of the age 
of meteorites from their helium content, and upon the gradual loss of 
the more volatile constituents of the meteors. The heating effect of 
the sun is increased through the character both of the systematic and 
accidental variation of g. A meteor met by the earth must have q 
less than 1, and because of the systematic increase of g with time the 
probability is that in the past history of the meteor the average q 
was still smaller. Furthermore the variation of g is such that small 
changes of g occur much more often, and are more symmetrically dis- 
tributed around the original value, than big perturbations. The 
result is a continuous unsteadiness of the orbit, a certain property of 
“mixing,” so that the probability is considerable that the meteor, 
before being caught by the earth, or before “jumping out” as the 
result of a big perturbation, will undergo some very close approaches 
to the sun. 

The stability of the orbit increases as the aphelion distance de- 
creases. At R = 2000 a. u., for a period of revolution of 30000 years, 
the perturbations during ¢ = 3.10° years are sensible, but not so 
revolutionary. Below 2000 a. u., the majority of the orbits may be 
regarded as statistically stable during the adopted interval of time. 


tity 
pig 
ee 
ak 
wif 
‘aap 
| 
4 


IS2 OPIK 


The perturbations of the inclination have much in common with 
the perturbations in perihelion distance except for the systematical 
character. When the perturbations are small, the average change of 


inclination is given by tgai = lait. for stronger perturbations the 


conditions are more complicated. In the following table is given the 
average expected change in inclination for t = 3.10° years. 


R 100 250 500 1000 1500 2000 
ai 414° +28° +55° + 82° + 110° 


For other intervals of time Ai is proportional to V t. An average 


change Ai = + 90° already means a close approach to uniform dis- 
tribution. According to statistics by Jantzen, 395 comets (1909) 
showed a more or less random distribution of the inclinations, with 
a slight preponderance of direct motions over retrograde ones; with 
an age of from 10° to 3.10% years, this would correspond to an average 
aphelion distance of 1500-2000 a. u., or a period of revolution of 
20000-30000 years. For greater aphelion distances the distribution 
of inclinations should be practically uniform, being smoothed out by 
perturbations. Meteors seem to show a much stronger preference for 
direct motion,! which could be preserved in the course of time only at 
a smaller effective aphelion distance, about 1000 a. u., corresponding 
to a period of about 10000 years. 


CONCLUSIONS. 


1. The possible limits of the solar system, as set by stellar pertur- 
bations during the time elapsed since the origin of the system, may be 
estimated at 10° astronomical units; whether the actual limits reach 
so far, depends upon the presence of bodies at these distances. 

2. Meteors, or comets of an aphelion distance exceeding 10000 
astronomical units, are not very likely to occur among the observable 
objects, because of the rapid increase of the average perihelion dis- 
tance due to stellar perturbations. 

3. From the scarce data relating to the observed nonuniformity of 
the orbital inclinations, the effective aphelion distance of observable 
comets is estimated at 1500-2000 astronomical units, and of observable 
sporadic meteors at 1000 astronomical units. 


1 E. Opik, H. C. 355, p. 9, 1930. 
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4. At the aphelion distance of 1000 a. u., the orbit of a meteor during 
3.10° years is comparatively stable; the average heating effect of the 
sun on an observable meteor corresponds in this case to from 100000 
to 300000 approaches at the actual perihelion distance, or to about 
30000 years of continuous heating to a temperature exceeding 0°C 
in vacuo, which may be sufficient for a stony meteorite to lose all its 
helium content. It is therefore probable that stony meteorites of 
solar origin yield very low ages according to the helium-radium 
method. 

5. As the result of stellar perturbations, the number of meteors and 
comets of a given constant aphelion distance of 2000 a. u. or greater 
must increase with increasing perihelion distance, to a certain limit; 
the rate of relative increase in the number must be greater for greater 
aphelion distances. For objects with an aphelion distance below 
1000 a. u., stellar perturbations could not change appreciably the 
original distribution of the orbital elements within 3.10° years. 

The present paper gives only an approximate sketch of the situa- 
tion. The approximate results are, however, quite sufficient to 
serve as bases for certain numerical estimates which appear to be 
useful and even necessary in the course of the systematic meteor 
research undertaken at Harvard Observatory under the initiative and 
organization of Dr. Shapley. 


Marcu 22, 1932. 
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